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Abstract. A simple re-scaling of velocities calculated from the Paczyhski & Wiita (1980) 
pseudo-Newtonian potential makes them consistent with special relativity and greatly im- 
proves the agreement with the exact relativistic calculations. The improvement is relevant 
in calculations of the Doppler effect, spectra and radiation transfer which all involve effects 
q of special relativity. 
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1. Introduction 

The pseudo-Newtonian potential (PN-potential) of Paczynski & Wiita (1980) provides 
a simple accurate model for black hole hydrodynamics. Under this approximation, one 
assumes that the gravitational potential of a non-rotating black hole with mass M is 

CM 2GM 

and one uses Newtonian theory. Here Rq is the gravitational radius of the black hole. The 
PN-potential has been used in many papers on black hole hydrodynamics, particularly on 
accretion disc dynamics (e.g., Okazaki, Kato, & Fukue 1987; Abramowicz et al. 1988). In 
short, the PN-potential is useful because it is both very simple and very accurate. For 
example, it predicts the same behaviour of the Keplerian circular orbits of free particles as 
the exact relativistic theory: around a non-rotating (Schwarzschild) black hole the orbits 
with R < 3Rq are unstable, and the orbits with R < 2Rq are unbound. Thus, the radius 
of the marginally stable orbit is located at R ms = 3Rg, and the radius of the marginally 
bound orbit at R m b = 2Rg- Accretion discs have their inner edges located somewhere 
between R ms and R m b- Conditions at radii smaller than R m b are not important for the 
astrophysical properties of accretion discs, but any Newtonian model of the black hole 
should be accurate for R > R m b- 

There is a problem with the approach based on the PN-potential. Even though the 
PN-potential describes accurately the dynamical effects of the curvature of spacetime, it 
does not include the effects of the special relativity. In particular, it does not guarantee 
that the velocity of matter should always be less than the speed of light. To see the 
problem, one may calculate the orbital Keplerian velocity of a test particle, v, and the 



corresponding Lorentz gamma factor, 7 = 1/^1 — (u/c) 2 , at the radius of the marginally 
bound orbit. The correct relativistic value is vbh = c/y/2 and the corresponding Lorentz 
gamma factor is jbh = \/2- However, the values calculated from the PN-potential are 
vpn = c and 7pjv = oo. Thus, it is quite obvious that the velocities and the Lorentz 
gamma factors are not well represented by the PN-potential. Since the Lorentz gamma 
factors are required in calculations of the Doppler effect, radiative transfer, and spectra, 
an improvement in calculating the Lorentz gamma factor is helpful in practice. 







2. Solution 

We suggest a re-scaling of the velocities calculated from the PN-potential which provides a 
simple solution to the problem. The velocity upjv calculated directly from the PN-potential 
should be replaced by Cpjv which satisfies 

Vpn = VpnJpn, 7pjv = , ■ (2) 

- (vpn/c) 2 

There is a "physical reason" for this re-scaling: if one writes relativistic formulae for 
the inertial forces in such a way that they resemble their Newtonian counterparts as much 
as possible, then v in the Newtonian formulae always corresponds to in the relativistic 
formulae (Abramowicz 1992). The re-scaling (2) follows from this general rule. More 
importantly, however, one may simply check that the re-scaling (2) works amazingly well in 
practice. In Figure la and b we show the velocities and the corresponding gamma factors 
for Keplerian circular orbits calculated according to the exact relativistic formulae, the 
PN-potential, and the PN-potential with the re-scaling (2). These velocities are expressed 



as 



PBH = = —====, (6 
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/3pjv = = \~ -, (4) 

c V z r — 1 



^ N = V = V2(r- r i)Hr' (5) 

where r = R/Rq. It is seen that, at r = 2, /3pjv = Pbh = The agreement between 

the re-scaled velocities and the exact relativistic ones is excellent for r > 2 as is seen from 
Figure lc and d. Specifically, for r > 2, vpn/vbh — 1 < 5% and 7pjv f^fBH — 1 < 2%. 

The same re-scaling may be used for other Newtonian models of black holes proposed 
more recently by several authors. All the other models are, however, less accurate and 
more complicated than that of Paczyhski and Wiita, so there is no obvious reason why 
one should use them (see Artemova, Bjornsson, & Novikov 1996) 

We conclude that the simple re-scaling of the velocities given by formula (2) makes the 
pseudo-Newtonian potential of Paczyhski and Wiita consistent with the special relativistic 
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effects. Agreement with the exact relativistic theory is better than ~ 5 percent. The re- 
scaled velocities calculated from the potential should be used in calculations of all special 
relativistic effects such as the relativistic Doppler shift, etc. 
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Figure 1. — The Keplerian orbital velocities (panel a) and the corresponding Lorentz 
gamma factors (panel b) calculated according to the exact relativistic formulae (solid line), 
the PN-potential (dotted line), and the PN-potential with the re-scaling (dashed line). The 
relative errors of the velocities (panel c) and the Lorentz gamma factors (panel d) calculated 
from the PN-potential with the re-scaling (dashed line) and without the re-scaling (dotted 
line). Note that at R = 2Rq } the PN-potential with the re-scaling gives the same result 
as that of the exact relativistic theory, and the agreement between these two is excellent 
for R > 2R G . 



